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Abstract 

Gradient-echo echo-planar imaging (GE EPI) is the most commonly used approach to assess localized blood oxygen level 
dependent (BOLD) signal changes in real-time. Alternatively, real-time spin-echo single-voxel spectroscopy (SE SVS) has 
recently been introduced for spatially specific BOLD neurofeedback at 3 T and at 7 T. However, currently it is not known 
how neurofeedback based on real-time SE SVS compares to real-time GE EPI-based. We therefore compared both methods 
at high (3 T) and at ultra-high (7 T) magnetic field strengths. We evaluated standard quality measures of both methods for 
signals originating from the motor cortex, the visual cortex, and for a neurofeedback condition. At 3 T, the data quality of 
the real-time SE SVS and GE EPI R2* estimates were comparable. At 7 T, the data quality of the real-time GE EPI acquisitions 
was superior compared to those of the real-time SE SVS. Despite the somehow lower data quality of real-time SE SVS 
compared to GE EPI at 7 T, SE SVS acquisitions might still be an interesting alternative. Real-time SE SVS allows for a direct 
and subject-specific T2* estimation and thus for a physiologically more plausible neurofeedback signal. 
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Introduction 

Gradient-echo echo-planar imaging (GE EPI) is the predomi- 
nant approach to assess localized blood oxygen level dependent 
(BOLD) signal changes [1]. Recent technological advances in the 
field of functional magnetic resonance imaging (fMRI) have made 
it possible to adapt GE EPI for use in real-time. Real-time fMRI 
can be used to provide a spatially localized neurofeedback signal 
based on the activation level of a region of the interest (ROI) [2— 
4], real-time brain-state classification [5-8], and connectivity- 
based neurofeedback [9]. Several studies have used real-time 
fMRI neurofeedback to train voluntary control over functionally 
specific brain areas, such as the motor and somatosensory cortices 
[10-13], the visual cortex [14,15], the cingulate cortex [8,16-19], 
the insula [20-22], the right inferior frontal gyrus [23], the 
amygdala [24-26], and the auditory cortex [27,28]. Some studies 
even suggested potential therapeutic effects of real-time fMRI 
neurofeedback training in chronic pain disorders [17], Parkinson's 
disease [29], tinnitus [28], schizophrenia [21,30], and depression 
[31]. 

Depending on the magnetic field strength, the acquisition 
technique, and the echo time, the intra- and the extravascular 
components contribute differently to the functional BOLD signal 
changes. Generally, the use of ultra-high magnetic field strengths 
increases the signal-to-noise ratio (SNR) of GE EPIs, as well as the 



sensitivity and specificity of the T2* BOLD contrast [32,33]. 
Higher magnetic field strengths also increase the SNR of the spin- 
echo (SE) EPI acquisitions which leads to improved localization of 
neural activity by targeting specifically microvasculature contri- 
bution to the BOLD signal [34,35]. This is typically achieved by 
SE EPI acquisitions and T2 contrast at high spatial resolution, i.e. 
by suppressing extra- and intravascular components from large 
vessels that normally contribute to the BOLD signal at GE [34,36- 
38] . Another difference between SE and GE techniques is that the 
maximum amplitude of the SE BOLD contrast is reached more 
quickly than that of the GE contrast. However, SE acquisitions 
usually have a lower SNR and a lower contrast-to-noise ratio 
(CNR) [35,39,40]. A disadvantage of ultra-high fields is that they 
are more prone to local field inhomogeneity, and thus require 
additional shimming adjustments and post-processing [41]. 

Functional SE single-voxel spectroscopy (SE SVS) has previ- 
ously been employed for neurofeedback applications by using a 
point-resolved spectroscopy (PRESS) acquisition protocol to 
acquire the large water peak in the spectrum [42—44]. Compared 
to conventional metabolite quantification methods where the 
water peak is typically suppressed [45-49], SE SVS acquisitions 
allow for a direct estimation of T2* from the unsuppressed water 
spectrum [44,50,51]. More specifically, the T2* neurofeedback 
signal was estimated directiy using free induction decay optimized 
linear regression, or using water peak non-linear parameterization 
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[42,43]. Such direct T2* estimates allow for localized and 
individually specific neurofeedback training, which is more 
physiologically plausible, and which extends beyond conventional 
neurofeedback techniques without direct T2* approximation. 
Another advantage of SE SVS acquisitions is that it allows for a 
reduction of specific absorption rate (SAR) levels, especially, as 
compared to SE EPI. This is particularly relevant for neurofeed- 
back training studies, where individuals are sometimes being 
scanned repeatedly for many hours. The SAR also significandy 
limits the number of slices that can be acquired with SE EPI 
sequences, especially at ultra-high magnetic fields [35]. Localized 
SE SVS is not restricted by these limitations once the neurofeed- 
back ROI has been defined. 

Overall, SE SVS acquisitions might be beneficial for neurofeed- 
back studies in that they allow for a direct and subject-specific 
T2*-based feedback signal [42], and in that they operate at lower 
SAR levels as compared to GE and SE EPI techniques. The 
possibility to simultaneously estimate T2* and T2 contrast [42] as 
well as reduced sensitivity to susceptibility-related signal loss of the 
SE techniques [35] provide an additional motivation to explore 
the SE SVS approach for neurofeedback. To shed light on these 
potential advantages, we for the first time provide a direct 
comparison analysis of R2* estimations based on GE EPI and SE 
SVS techniques at 3 T and at 7 T magnetic fields. Our analysis 
involved standard quality measures such as percent signal change 
(A%), contrast-to-noise ratio (CNR), and t-statistics for region- 
specific time courses during standard functional localizer runs and 
during neurofeedback runs. Notably, we did not compare signal 
quality between 3 T and 7 T magnetic field strengths, which have 
previously been thoroughly investigated for SE EPI as well as for 
GE EPI techniques. 

Methods 

Experimental Design 

Two different groups of 7 healthy volunteers each were scanned 
on a 3 T scanner (4 male, 3 female, age 28 ±7 years) and on a 7 T 
scanner (6 male, 1 female, age 33±9 years). All participants were 
right-handed according to a minimal score of 6 on the Edinburgh 
Handedness Inventory [52] . Study protocols were approved by the 
Ethics Committees of the Medical Faculty of the RWTH Aachen 
University and of the University of Pennsylvania. All participants 
gave written informed consent and were paid an allowance at the 
end of their participation. 

The experimental protocol which was used to compare real- 
time GE EPI and real-time SE SVS acquisitions consisted of the 6 
following runs (Figure 1): a GE EPI and a SE SVS based 
functional localizer of the primary motor cortex (PMC loc), a SE 
SVS and a GE EPI based neurofeedback run targeting the 
primary motor cortex (PMC NF), and a GE EPI and a SE SVS 
based localizer of the visual cortex (VC loc). 

The primary motor cortex functional localizer runs consisted of 
finger tapping and baseline blocks, and the visual cortex functional 
localizer runs consisted of the presentation of a flickering visual 
checkerboard and baseline blocks. For the neurofeedback runs, the 
participants were instructed to adjust the speed and strength of 



their finger tapping so that a green horizontal bar would move up 
to the level of a predefined red horizontal target bar. All functional 
runs comprised 5 blocks of activation (i.e. finger tapping or visual 
stimulation, respectively), interleaved with 5 baseline blocks. Each 
block lasted 30 seconds, resulting in total run duration of 5 min. 

Data Acquisition 

Functional GE EPI and SE SVS data were acquired on a 3 T 
and a 7 T MR scanner (Siemens Medical Solutions, Erlangen, 
Germany) equipped with a transmit body coil and a 12-channel 
phased array head receive coil at 3 T, or a birdcage single-channel 
head coil (quadrature) at 7 T. EPI images were obtained with a 
single-shot gradient-echo T2*-weighted sequence with 300 repe- 
titions (TR= 1000 ms, 16 slices, volumes matrix size 64x64, voxel 
size = 3 x3 x3.75 mm 3 , flip angle ot = 77°, bandwidth = 2.23 kHz/ 
pixel, TE = 30 ms at 3 T; TE = 28 ms at 7 T). At 3 T, the water 
spectra were acquired using a spin-echo PRESS protocol with 300 
repetitions (TE/TR= 30/1000 ms, flip angle ot = 90°- 180°- 180°, 
bandwidth = 1 kHz, acquisition duration = 512 ms). At 7 T, the 
acquisition protocol was slightly different with TE = 20 ms, 
bandwidth = 2 kHz, acquisition duration = 256 ms. Spectroscopic 
voxels were chosen as isotropic as possible based on the individual 
GE EPI brain activation maps for the motor and visual conditions 
(approximately lxlxl cm' 1 ). On both scanners we performed a 
manual calibration of the transmitter amplitude and optimized the 
gradient shim currents using Siemens manual shimming adjust- 
ments in order to improve the spectroscopic signal quality. 
Acquisition parameters were selected to obtain a robust T2* 
estimate and the BOLD effect. Note that for the SE SVS pulse 
sequence, the TE was selected as short as possible for the given 
PRESS protocols, i.e. we targeted the T2* contrast. This was also 
done to reduce T2 weighting, and to acquire early-echo data for a 
more accurate T2* approximation. The SE SVS estimates of the 
T2* contrast were barely affected by the flip angle for the 
transversal magnetization. This was because the inversion pulses 
contributed to the Tl saturation, because the post-acquisition 
delay time was long compared to the T2 of the tissue, and because 
the T2* estimates were calculated from the free induction decay 
function (FID) direcdy. Neither water nor fat suppression was 
applied for the spectroscopy protocols. The first 10 acquisitions 
were discarded to avoid Tl saturation effects. The visual 
instructions and feedback were shown to the subjects via MR- 
compatible goggles (Resonance Technology Inc. Northridge USA) 
on the 3 T scanner, and projected to an MR-compatible screen on 
the 7 T scanner. The data were exported in real-time to the local 
PC and processed with the custom-made software as described in 
[43,53]. 

Data Processing and Feedback Signal Extraction 

Immediately after acquiring the data from the primary motor 
and visual cortex localizer runs, the images were pre-processed 
with SPM8 functions (Wellcome Trust Centre for Neuroimaging, 
Queen Square, London, UK), i.e. realigned to the first scan of the 
respective localizer run, and smoothed with an isotropic Gaussian 
kernel with 4 mm full-width-at-half-maximum (FWHM). Next, we 
specified a general linear model (GLM) with regressors for the 
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Figure 1. Sequence of data acquisition. GE - gradient-echo, EPI - echo-planar imaging, SE - spin-echo, SVS - single voxel spectroscopy, PMC - 
primary motor cortex, VC - visual cortex, NF - neurofeedback, loc - localizer. 
doi:1 0.1 371 /journal.pone.0091 620.g001 
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Figure 2. Illustration of the GE EPI and SE SVS ROIs. A GE EPI activation map and a single-voxel PMC ROI (blue) of a representative participant 
are shown on sagittal, transverse, and coronal planes of this participant's structural scan. The PMC SE SVS ROI (size approximately 1x1x1 cm 3 ) was 
defined to cover the voxels that exhibited a significant positive BOLD response to the left hand finger tapping. The GE EPI ROI was restricted to the SE 
SVS ROI. 

doi:1 0.1 371 /journal.pone.0091 620.g002 

experimental conditions, and acquired thresholded t-maps. For response to the left hand finger tapping or visual stimulation, 
the SE SVS acquisitions, the single-voxel was defined so that it respectively [43]. The ROI for the GE EPI acquisitions was 
covered those voxels that exhibited a significant positive BOLD restricted to the SE SVS ROI (Figure 2). 




time, s time, s 

Figure 3. T2* approximation of the SE SVS data acquired at 3 T and at 7 T. The linear regression fits (red) are shown for single PMC ln|FID|'s 
(blue) for representative participants at 3 T (A) and at 7 T (B). Optimal linearization lengths are 0.18 s at 3 T (t = 41.6, p<0.001), and 0.052 s at 7 T 
(t = 1 6.7, p<0.001). 

doi:1 0.1 371 /journal.pone.0091 620.g003 
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Figure 4. Block-related averages of the GE EPI and the SE SVS R2* time courses. Block-related averages were calculated for the 3 T GE EPI 

(blue), the 3 T SE SVS (green), the 7 T GE EPI (black), and the 7 T SE SVS (red) time courses of the PMC, the PMC NF, and the VC runs. Error bars 
represent the standard error of the mean. 
doi:1 0.1 371 /journal.pone.0091 620.g004 



During the neurofeedback primary motor cortex GE EPI runs, 
the EPI volumes were first realigned to the first volume of the 
primary motor cortex functional localizer run. The feedback 
signal, which corresponded to the average activity within the ROI, 
was then calculated as soon as a new volume was acquired. During 
the neurofeedback SE SVS runs, the acquired water spectra were 
shifted to zero, filtered with a Gaussian filter, the eddy currents 
were compensated, and the water spectra were phase-corrected 
[54]. The feedback was provided after each FID acquisition as an 
absolute T2* measure which was estimated with the statistically 
optimized linear regression approach [43]. The optimal linear 
regression length was estimated based on the SE SVS primary 
motor and visual cortex runs. After the feedback signal was 
extracted from either SE SVS or GE EPI acquisitions, the signal 
was processed in order to reduce noise and to remove spike-like 
artifacts using our custom-made real-time software [53]. For the 
GE EPI acquisitions, the head motion parameters were taken into 
account, but head motion parameters were not available for the 
SE SVS acquisitions. To ensure that motion artifacts did not cause 
significant SE SVS signal distortions, we located relatively small 
ROIs within large active zones revealed by the primary motor and 
visual cortex localizer runs. Inter-run head movements between 
the SE SVS runs were controlled by acquiring GE EPI scans 
before and after the SE SVS runs; they were less than 1 mm. 

Time Courses Quality Measures and Comparison Analysis 

The comparison analysis between GE EPI and SE SVS 
acquisitions was based on their CNR, percent signal change, 
and t-statistics, and was performed separately for data acquired at 
3 T and at 7 T. For SE SVS, statistically optimized linear 
regression was applied to the acquired FID in the time domain. 
The natural logarithm of the FID can be simplified assuming that 
the water signal is the dominating component in the acquired FID, 
and that all other proton sources of the signal are negligible [43]: 



In (FID) = - t/T2* water + In (A water ) 



curve {ln(\FID\)) and determined by the optimal linear regression 
length (OLR) [43]. To compensate for line-broadening caused by 
applied Gaussian nitration, T2* estimation function was weighted 
with the correspondent filter coefficients [42]. The optimal linear 
regression length was estimated in the sense of a statistical 
measure, i.e. the maximum t-value in the distribution of t-values of 
time series computed for a set of linear regression lengths (Figure 3; 
red curves). The processed signal in Equation [1] may still have a 
large non-linear component because of the inadequate shimming 
conditions (Figure 3; blue curves), which can complicate the 
regression analysis of the acquired data. However, despite its 
simplicity, the proposed OLR approach has been shown to 
provide reliable T2* estimations at high and ultra-high magnetic 
fields [42,43]. 

Because the GE EPI voxel intensity is proportional to exp(— 
TE-R2*), the T2* values estimated from SE SVS time courses were 
transformed to exp( — TE/T2*) using the applied echo time TE 
and arbitrary scaling. This allowed for a direct comparison 
between GE EPI and SE SVS acquisitions. 

Block-related averages were averaged across the time-course 
condition/baseline periods. The percent signal changes (A%) were 
estimated as an average from block-related condition/baseline 30- 
point plateaus. For the statistical analysis of the BOLD signal 
changes, we specified general linear models (GLM) with regressors 
for the experimental conditions defined in SPM8 (W elcome Trust 
Centre for Neuroimaging, UK). Each participant's fMRI motion 
parameters were included into the GLM as nuisance regressors. 
Effects on the time-course quality ratings were analyzed in 
repeated-measures ANOVA for all data sets with functional run 
(PMC, PMC NF, and VC), MR scanner (3 T and 7 T) and 
acquisition technique (GE EPI and SE SVS) as within-subject 
factors. To further evaluate the difference between two samples, 
standard two-sample t-tests were used (t- and p- values; one-tailed). 
To calculate the CNR, we estimated differences between signal 
means during baseline and activation blocks, and their residual 
variances: 



with water time constant T2* and amplitude A. The linear 
regression was subsequently applied to the absolute logarithmic 



CNR-- 



mean(condition) — mean(baseline) 
\J \ax(condition) + \ax(baseline) 



(2) 
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Figure 5. Performance comparison between GE EPI and SE SVS time courses at 3 T and at 7 T. We compared the time course quality of 
the PMC (1 st column), the PMC NF (2 nd column), and the VC (3 rd column) in terms of their percent signal change (A%; panel A), contrast-to-noise ratio 
(CNR, panel B), and t-statistics (t-value, panel C). This was done separately for data acquired with GE EPI (blue bars) and SE SVS (red bars), and 
separately for 3 T and 7 T acquisitions. At 3 T, the higher SE SVS CNR and t-value were indicated with red arrows. At 7 T, the higher GE EPI A% was 
indicated with green arrows. Error bars represent the standard error of the mean; asterisks denote statistical significance (p<0.05). 
doi:1 0.1 371 /journal.pone.0091 620.g005 
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where condition /baseline is the time course of the ROI in the 
functional localizer condition and baseline, respectively. All 
computations were carried out on a standard PC in Matlab 7.10 
(The Mathworks, Natick, MA). The custom-made neurofeedback 
toolbox is available on request from the corresponding author. 

Results 

To illustrate the quality of the acquired time courses, we 
evaluated the BOLD-dependent block-related signal changes in 
the GE EPI and the SE SVS time courses (Figure 4). 

Overall, GE EPI and SE SVS acquisitions showed high data 
quality at 3 T and at 7 T in terms of the applied quality measures 
(Figure 5, see also Table 1 for numeric values). Compared to our 
previous study [42] , where SE SVS time courses were estimated in 
terms of the T2* measures, R2*-weighting of the time-courses led 
to similar results. We found that the average R2* values in the 
PMC (18.9±0.1 s" 1 ) and in the PMC NF runs (18.5±0.2 s" 1 ) at 
3 T, in the VC runs at 3 T (24.6±0.3 s" 1 ), and in the VC runs at 
7 T (55.9±0.5 s ) were similar to previous findings [38,32,55]. 
However, the R2* values in the PMC (61.0±0.7 s" 1 ) and in the 
PMC NF (59.8±0.9 s" 1 ) runs at 7 T were somewhat higher. A 
repeated-measures ANOVA revealed a significant interaction of 
the factors scanner*technique (A%: F(l,12) = 5.7, p = 0.034; CNR: 
F(l,12) = 4.90, p = 0.047; t-statistics: F(l,12) = 5.7, p = 0.034). In 
addition, for percentage signal change, the interaction of the 
factors scanner*run was significant (F(2,24) = 4. 1, p = 0.029). This 
implied that neither the data acquisition technique nor the field 
strength appeared to have an unambiguous advantage. Instead the 
performance depended on the specific combination of field 
strength, acquisition technique, and ROI. We therefore evaluated 
these factors using pair-wise comparisons. 

At 3 T, the GE EPI and the SE SVS acquisitions did not differ 
significandy in terms of A% for any of the runs (Figure 5 A, PMC, 
PMC NF, VC: all t<0.8, all p>0.2). In contrast, at 7 T, A% was 
significandy higher for GE EPI acquisitions as compared to SE 
SVS acquisitions in the PMC (t= 1.8, p = 0.046), and in the VC 
runs (t = 1 .9, p = 0.043). The same trend, albeit non-significant was 
evident also for the PMC NF runs (t = 0.8, p>0.2). 

Comparing GE EPI and SE SVS acquisitions at 3 T in terms of 
their CNR (Figure 5B) showed significantly higher differences in 
the PMC runs (t = 2.1, p = 0.027), and smaller differences in the 
PMC NF (t = 1 . 1 , p = 0. 1 5) and in the VC runs (t = 1 .5, p = 0.079). 
At 7 T, there was no significant CNR difference between the GE 
EPI and SE SVS acquisitions (PMC, PMC NF, VC: aU t <0.8, all 
p>0.2). Notably, the SE SVS acquisitions at 7 T showed larger 
variations in the applied quantitative measures (Table 1). 

The pattern of results for the t-statistics was similar to those of 
the CNR (Figure 5C), which is not surprising because they 
represent similar metrics. At 3 T, the t-values of the SE SVS 
acquisitions were significantly higher than those of the GE EPI 
acquisitions for the PMC NF runs (t = 2.7; p = 0.01), and higher 
for the PMC runs (t= 1.6; p = 0.067). No other differences 
between GE EPI and SE SVS acquisitions were found. 

Discussion 

GE EPI vs. SE SVS Analysis 

We showed that real-time SE SVS at 3 T led to comparable 
data quality with that of real-time GE EPI. The BOLD percent 
signal changes for the GE EPI and for the SE SVS acquisitions at 
3 T were stable, physiologically plausible, and did not differ 
between the acquisition methods (Figures 4, 5). Also, the CNR and 
the t-statistics of the GE EPI and the SE SVS acquisitions were 
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comparable at 3 T. In some conditions, the SE SVS acquisitions at 
3 T even showed enhanced performance compared to that of the 
GE EPI acquisitions considering the voxel size usually applied in 
real-time fMRI studies (Figure 5; red arrows). 

At 7 T, the BOLD percent signal changes increased for both 
acquisition methods, but the signal changes of the GE EPI 
acquisitions were higher than that observed for the SE SVS 
acquisitions. Although our acquisition protocol was optimized for 
the T2* contrast, the estimation at 7 T was affected by the T2 
contrast even at the shortest TE possible for the given MR 
sequence [39]. 

On the other hand, real-time SE SVS might allow for providing 
the T2 contrast at the same time as the T2* contrast in order to 
target specifically the microvasculature [37,38,42]. The decay rate 
of the FID is weighted by the intensity at time zero (To), which is 
followed by the T2/T1 contrast: \FID\=I 0 exp(-t/T2*). Note, 
that in our SE SVS pulse sequence, the SE data acquisition starts 
at time TE (i.e. time zero) and the percent T2 changes can be 
approximated while neglecting the Tl effect, i.e. 
T2~ — TE/\n(Io) if Iq is given. Taking into account that the 
latter estimation was applied for a shorter than canonical spin- 
echo TE used for T2 contrast, and that it could be biased if the 
FID contains multiple components, the calculated percent signal 
changes values were ~0.2% [42]. Note, that the present study was 
designed for an optimal T2* approximation and could be further 
balanced for higher T2 contrast, e.g. by using a larger TE 
[37,38,56,57]. This supports that real-time SE SVS with longer 
TEs might allow for a neurofeedback signal originating from 
specifically the microvasculature T2 and T2* dynamics. Further 
research is needed to explore this possibility. 

The A% estimated at 7 T confirmed that GE EPI yields a very 
high BOLD contrast at ultra-high field strengths. However, the 
advantage of ultra-high magnetic fields is reflected in improved 
CNR and improved statistics only if noise level can be controlled; 
particularly increased susceptibility artifacts, and local field 
inhomogeneity. In our study, the signal change increased at 7 T, 
but not the CNR and the t-statistics, which indicates that noise 
level increased as well. This is also illustrated by the increased 
variance of the R2* time courses (Figure 4). Also, whereas the 
average R2* values in the VC runs at 3 T and at 7 T, and in the 
PMC runs at 3 T were plausible, the R2* values in the PMC runs 
at 7 T were somewhat higher. The latter and the fact that the T2* 
approximations were more stable in the VC runs [42,43] suggest 
suboptimal shimming in the PMC runs. Recendy proposed local 
Bi (Bi" 1 ") shimming in combination with a multichannel transceiver 
array coil might address this limitation [58]. Additionally, static 
magnetic field inhomogeneity (B 0 ) for a single-voxel approach can 
be relatively easy compensated with a strong second-order and, if 
available, third-order shim system [47]. 

In our study, the difference between acquisitions did not only 
depend on the acquisition method, but was also region specific. 
For example, reduced CNR and t-statistics for the GE EPI runs at 
7 T was also observed for the VC ROI, but not for the PMC ROI 
(Figure 5, Table 1). This might be due to the suboptimal shimming 
and, for SE SVS, the linearization length function (i.e. the 
individual T2*/R2* estimation), which can be very specific 
depending on the ROI and on the magnetic field strength 
[42,43]. Interestingly, our GE EPI PMC and PMC NF runs 
benefit more from 7 T compared to the VC runs (Figure 5). 
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acquisition might allow for providing a combined T2* and T2 
contrast, and thus potentially for targeting specifically the macro- 
and microvasculature, for reducing SAR for scanning-intensive 
neurofeedback training, and for reducing sensitivity to suscepti- 
bility-related signal loss in specific brain regions. However, these 
potential advantages need to be experimentally validated in future 
studies. 

Author Contributions 

Conceived and designed the experiments: YK ME KM. Performed the 
experiments: YK. Analyzed the data: YK ME KM. Contributed reagents/ 
materials/analysis tools: YK ME FS KM. Wrote the paper: YK ME FS 
KM. 



2. dcCharms RC (2008) Applications of real-time fMRI. Nat Rev Neurosci 9: 720- 
729. 



PLOS ONE | www.plosone.org 



7 



March 2014 | Volume 9 | Issue 3 | e91620 



Real-Time Water Proton SE PRESS vs GE EPI at 3 & 7 T 



3. Weiskopf N (2012) Real-time 1MRI and its application to neurofeedback. 
Neuroimage 62: 682-692. 

4. SulzerJ, Haller S, Schamowski F, Weiskopf N, Birbaumcr N, ct al. (2013) Real- 
time fMRI neurofeedback: Progress and challenges. Neuroimage 76: 386-399. 

5. LaContc SM, Peltier SJ & Hu XP (2007) Real-time fMRI using brain-state 
classification. Hum Brain Mapp 28: 1033-1044. 

6. LaContc SM (201 1) Decoding 1MRI brain states in real-time. Neuroimage 56: 
440-454. 

7. Sitaram R, Weiskopf N, Caria A, Veit R, Erb M, et al. (2008) 1MRI brain- 
computer interfaces. IEEE Signal processing magazine 25: 95—106. 

8. Sitaram R, Lee S, Ruiz S, Rana M, Veit R, ct al. (2011) Real-time support 
vector classification and feedback of multiple emotional brain states. Neuro- 
image 56: 753-765. 

9. Koush Y, Rosa MJ, Robincau F, Hcincn K, Ricgcr S, et al. (2013) Connectivity- 
based neurofeedback: dynamic causal modeling for real-time fMRI. Neuroimage 
81: 422-430. 

10. Chicw M, LaContc SM & Graham SJ (2012) Investigation of fMRI neurofeed- 
back of differential primary motor cortex activity using kinesthetic motor 
imagery. Neuroimage 61: 21—31. 

11. Sitaram R, Veit R, Stevens B, Caria A, Gerloff C, et al. (2012) Acquired control 
of ventral premotor cortex activity by feedback training: an exploratory real-time 
FMRI and TMS study. Neurorehabil Neural Repair 26: 256-265. 

12. dcCharms RC, Christoff K, Glover GH, PaulyJM, Whitfield S, et al. (2004) 
Learned regulation of spatially localized brain activation using real-time fMRI. 
Neuroimage 21: 436-443. 

13. Bray S, Shimojo JP, G'Doherty JP (2007) Direct instrumental conditioning of 
neural activity using functional magnetic resonance imaging-derived reward 
feedback. J Ncurosci 27: 7498-7507. 

14. Scharnowski F, Hutton C, Josephs G, Weiskopf N, Recs G (2012) Improving 
visual perception through neurofeedback. J Ncurosci 32: 17830-17841. 

15. Shibata K, Watanabe T, Sasaki Y, Kawato M (2011) Perceptual Learning 
Incepted by Decoded fMRI Neurofeedback Without Stimulus Presentation. 
Science 9: 1413-1415. 

16. Weiskopf N, Veit R, Erb M, Mathiak K, Grodd W, ct al. (2003) Physiological 
self-regulation of regional brain activity using real-time functional magnetic 
resonance imaging (fMRI): methodology and exemplary data. Neuroimage 19: 
577-586. 

17. dcCharms RC, Macda F, Glover GH, Ludlow D, PaulyJM, ct al. (2005) 
Control over brain activation and pain learned by using real-time functional 
MRI. PNAS 102: 18626-18631. 

18. Mathiak KA, Koush Y, Dyck M, Gabcr TJ, Alawi E, ct al. (2010) Social 
reinforcement can regulate localized brain activity. Eur Arch Psychiatry Clin 
Ncurosci 260: 132-136. 

19. Hamilton JP, Glover GH, Hsu JJ, Johnson RF & Gotlib IH (201 1) Modulation 
of Subgenual Anterior Cingulate Cortex Activity With Real-Time Neurofeed- 
back. Human Brain Mapping 32: 22-31. 

20. Veit R, Singh V, Sitaram R, Caria A, Rauss K, ct al. (2012) Using real-time 
fMRI to learn voluntary regulation of the anterior insula in the presence of 
threat-related stimuli. Social Cognitive and Affective Neuroseicnce 7: 623—634. 

21. Ruiz SM, Sitaram R, Lee S, Sockadar S, Birbaumcr N (2010) Learning to self- 
regulate insula cortex modulates emotion recognition and neural connectivity in 
schizophrenia. Schizophrenia Research 117: 178. 

22. Caria A, Veit R, Sitaram R, Lotzc M, Weiskopf N, et al. (2007) Regulation of 
anterior insular cortex activity using real— time 1MRI. Neuroimage 35: 1238— 
1246. 

23. Rota G, Sitaram R, Veit R, Erb M, Weiskopf N, ct al. (2009) Self-regulation of 
regional cortical activity using real-time fMRI: the right inferior frontal gyrus 
and linguistic processing. Human Brain Mapping 30: 1605-1614. 

24. Posse S, Fitzgerald D, Gao KX, Habel U, Rosenberg D, et al. (2003) Real-time 
fMRI of temporolimbic regions detects amygdala activation during single-trial 
self-induced sadness. Neuroimage 18: 760-768. 

25. Johnston SJ, Boehm SG, Healy D, Goebel R, Linden DEJ (2010) Neurofeed- 
back: A promising tool for the self-regulation of emotion networks. Neuroimage 
49: 1066-1072. 

26. Zotcv V, Krucgcr F, Phillips R, Alvarez RP, Simmons WK, ct al. (2011) Self- 
Regulation of Amygdala Activation Using Real-Time fMRI Neurofeedback. 
PloS ONE 6: e24522. 

27. Yoo SS, LeeJH, O'Leary H, Lec V, Choo SE, ct al. (2007) Functional magnetic 
resonance imaging-mediated learning of increased activity in auditory areas. 
Ncurorcport 18: 1915-1920. 

28. Haller S, Birbaumcr N & Veit R (2010) Real-time fMRI feedback training may 
improve chronic tinnitus. European Radiology 20: 696-703. 

29. Subramanian L, Hindlc JV, Johnston S, Roberts MV, Husain M, ct al. (201 1) 
Real-time functional magnetic resonance imaging neurofeedback for treatment 
of Parkinson's disease. J Ncurosci 31: 16309-16317. 

30. Ruiz S, Lee S, Sockadar SR, Caria A, Veit R, et al. (2013) Acquired Self-control 
of Insula Cortex Modulates Emotion Recognition and Brain Network 
Connectivity in Schizophrenia. Hum Brain Mapp 34: 200-212. 

31. Linden DEJ, Habcs I, Johnston SJ, Linden S, Tatineni R, et al. (2012) Real- 
Time Sclf : Rcgulation of Emotion Networks in Patients with Depression. PloS 
ONE7:c38115. 

32. Van Dcr Zwaag W, Francis S, Head K, Peters A, Gowland P, ct al. (2009) fMRI 
at 1.5, 3 and 7 T: Characterising BOLD signal changes. Neuroimage 47: 1425- 
1434. 



33. Declchand DK, Van De Moortelc PV, Adriany G, litis I, Andersen P, et al. 
(2010) In vivo 1H NMR spectroscopy of the human brain at 9.4 T: Initial 
results. Journal of Magnetic Resonance 206: 74—80. 

34. Ugurbil K, Adriany G, Andersen P, Chen W, Garwood M, ct al. (2003) 
Ultrahigh field magnetic resonance imaging and spectroscopy. Magn Reson 
Imaging 21: 1263-1281. 

35. Budde J, Shajan G, Zaitsev M, Schefller K and Pohmann R (2013) Functional 
MRI in Human Subjects with Gradient-Echo and Spin-Echo EPI at 9.4 T. 
Magn Reson Med 71: 209-218. 

36. Yaeoub E, Schmucl A, Logothetis N, Ugurbil K (2007) Robust detection of 
ocular dominance columns in humans using Halm Spin Echo BOLD functional 
MRI at 7 Tesla. Neuroimage 4: 1161-1177. 

37. Yaeoub E, Duong TQ Van de Moortele PF, Lindquist M, Adriani G, et al. 
(2003) Spin-echo fMRI in humans using high spatial resolutions and high 
magnetic fields. Magn Reson Med 49: 655-664. 

38. Duong TQ, Yaeoub E, Adriany G, Hu X, Ugurbil K, ct al. (2003) 
Microvascular BOLD contribution at 4 and 7 T in the human brain: 
gradient-echo and spin-echo fMRI with suppression of blood effects. Magn 
Reson Med 49: 1019-1102. 

39. Hulvcrshorn J, Bloy L, Gualtieri EE, Leigh JS, and Elliott MA (2005) Spatial 
sensitivity and temporal response of spin echo and gradient echo bold contrast at 
3 T using peak hemodynamic activation time. Neuroimage 24: 216-223. 

40. Schmidt CF, Boesiger F Ishai A (2005) Comparison of fMRI activation as 
measured with gradient- and spin-echo EPI during visual perception, Neuro- 
image 26: 852-859. 

41. Tkac I and Gruetter R (2005) Methodology of H NMR Spectroscopy of the 
Human Brain at Very High Magnetic Fields. Appl Magn Reson 29(1): 139-157. 

42. Koush Y, Elliott MA, Scharnowski F, Mathiak K (2013) Real-time automated 
spectral assessment of the BOLD response for neurofeedback at 3 and 7 T. 
J Neurosci Methods 218: 148-160. 

43. Koush Y, Elliott MA, Mathiak K (201 1) Single Voxel Proton Spectroscopy for 
Neurofeedback at 7 Tesla. Materials 4: 1548-1563. 

44. Zhu X and Chen W (2001) Observed BOLD Effects on Cerebral Metabolite 
Resonances in Human Visual Cortex During Visual Stimulation: A Functional 
1H MRS Study at 4 T. Magn Reson Med 46: 841-847. 

45. Schaller BM, Mekle R, Xin L, O'brien K, Gruetter R (2013) Net Increase of 
Lactate and Glutamate Concentration in Activated Human Visual Cortex 
Detected With Magnetic Resonance Spectroscopy at 7 Tesla. J of Neurosci Res 
91: 1076-1083. 

46. Mangia S, Tkac I, Gruetter R, Van De Moortele PF, Giovc F, ct al. (2006) 
Sensitivity of single-voxel 1H-MRS in investigating the metabolism of the 
activated human visual cortex at 7 T. Magn Reson Imaging 24: 343-348. 

47. Tkac I, Oz G, Adriany G, Ugurbil K and Gruetter R (2009) In Vivo 1H NMR 
Spectroscopy of the Human Brain at High Magnetic Fields: Metabolite 
Quantification at 4T vs. 7T. Magn Reson Med 62: 868-879. 

48. Van dcr Veen JWC, Weinberger DR, Tcdcschi G, Frank JA, Duyn JH (2000) 
Proton MR Spectroscopic Imaging without Water Suppression. Radiology 217: 
296-300. 

49. Haasc A, Frahm J, Hacnickc W, Matthci D (1985) 1H NMR chemical shift 
selective (CHESS) imaging. Phys Med Biol 30: 341-344. 

50. Richards TL, Dager SR, Posse S (1998) Functional MR spectroscopy of the 
brain. Neuroimaging Clin N Am 8: 823-824. 

51. Hennig J, Ernst T, Speck O (1994) Detection of brain activation using 
oxygenation sensitive functional spectroscopy. Magn Reson Med 31: 85—90. 

52. Oldfield RC (1971) The assessment and analysis of handedness: the Edinburgh 
inventory. Neuropsychologia 9: 97—113. 

53. Koush Y, Zvyagintsev M, Dyck M, Mathiak KA, Mathiak K (2012) Signal 
quality and Bayesian signal processing in neurofeedback based on real-time 
(MRI, Neuroimage 59: 478-489. 

54. Klose U (1990) In Vivo Proton Spectroscopy in Presence of Eddy Currents, 
Magn Reson Med 14: 26-30. 

55. Peters AM, Brookes MJ, Hoogenraad FG, Gowland PA, Francis ST, ct al. (2007) 
T2* measurements in human brain at 1.5, 3 and 7 T. Magn Reson Imaging 25: 
748-753. 

56. Silvcnnoinen MJ, Clingman CS, Golay X, Kauppincn RA, van Zijl PC (2003) 
Comparison of the dependence of blood R2 and R2* on oxygen saturation at 1.5 
and 4.7 Tesla. Magn Reson Med 49: 47-60. 

57. Nair G, Duong TQ(2004) Echo-planar BOLD fMRI of mice on a narrow-bore 
9.4 T magnet. Magn Reon Med 52: 430-434. 

58. Emir UE, Auerbach EJ, Van De Moortele PF, Marjanska M, Ugurbil K, ct al. 
(2012) Regional neurochemical profiles in the human brain measured by 1H 
MRS at 7 T using local Bl shimming. NMR Biomed 25: 152-160. 

59. Alger JR (2010) Quantitative Proton Magnetic Resonance Spectroscopy and 
Spectroscopic Imaging of the Brain: A Didactic Review. Top Magn Reson 
Imaging 21: 115-128. 

60. Posse S, Otazo R, Tsai SY, Yoshimoto AE, Lin FH (2009) Single-Shot Magnetic 
Resonance Spectroscopic Imaging with Partial Parallel Imaging. Magn Reson 
Med 61: 541-547. 

61. Posse S, Otazo R, Caprihan A, Bustillo J, Chen H, ct al. (2007) Proton echo- 
planar spectroscopic imaging ofj-coupled resonances in human brain at 3 and 4 
Tesla. Magn Reson Med 58: 236-244. 

62. Sijens PE, Dorrius MD, Kappert P, Baron P, Pijnappcl RM, ct al. (2010) 
Quantitative multivoxel proton chemical shift imaging of the breast. Magn 
Reson Imaging 28: 314-319. 



PLOS ONE | www.plosone.org 



8 



March 2014 | Volume 9 | Issue 3 | e91620 



Real-Time Water Proton SE PRESS vs GE EPI at 3 & 7 T 



63. Mulkern RV, Chen N, Oshio K, Panych LP, Rybicki FJ, et al. (2004) Fast 
spectroscopic imaging strategics tor potential imaging applications in 1MRI. 
Magn Reson Imaging 22: 1395-1405. 

64. Tsai SY, Posse S, Lin YR, Ko CW, Otazo R, et al. (2007) Fast Mapping of the 
T2 Relaxation Time of Cerebral Metabolites Using Proton Echo-Planar 
Spectroscopic Imaging (PEPSI). Magn Reson Med 57: 859-865. 



65. Maudslcy AA, Hilal SK, Pcrman WH, Simon HE (1983) Spatially resolved high- 
resolution spectroscopy by 'four-dimensional' NMR. J Magn Reson 51: 147— 
152. 

66. Wciskopf N, Klosc U, Birbaumcr N, Mathiak K (2005) Single-shot 
compensation of image distortions and BOLD contrast optimization using 
multi-echo EPI for real-time fMRI. Ncurolmagc 24: 1068-1079. 



PLOS ONE | www.plosone.org 



9 



March 2014 | Volume 9 | Issue 3 | e91620 



